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Abstract 

In models with large extra dimensions, where quantum gravity effects be- 
come strong at the TeV scale, the p parameter can receive large contributions 
from one-loop diagrams involving exchange of multiple graviton and dilaton 
states. These contributions are computed, taking into account cancellation 
of spurious infrared divergences, and the (finite) results for d = 5 and 6 are 
compared with current experimental data. It is shown that 5 large extra 
dimensions are incompatible with the data and d = 6 is severely constrained. 
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1. Introduction 



Considerable excitement has been generated by the suggestion [|T| that strong gravita- 
tional effects could become manifest at scales of the order of a few TeV. One of the ways 
in which this can happen requires large compactified spacetime dimensions beyond the 
four Minkowski dimensions. An elegant model involves a string theory 0] (living in 10 di- 
mensions) which has solitonic excitations of the gravitational field called D-branes 0. In 
this model, proposed Q by Antoniadis, Arkani-Hamed, Dimopoulos and Dvali (AADD), 
all Standard Model (SM) fields are conceived of as living on a heavy D^-hrane embedded 
in a space of 10 dimensions. In fact, the SM fields correspond to excitations of open 
strings whose ends are confined to the brane. Of the extra 6 dimensions, it is possible to 
have d of them compactified with radii Rc, where Rc can be large, while the remaining 
(6-(i) dimensions are compactified with radii near the Planck length Mpl ~ 10"'^^ cm. 
The latter do not play much part in the subsequent discussion, though their existence is 
essential to have a consistent string theory in the first place. For all practical purposes, 
therefore, spacetime consists of (4 + d) dimensions, the extra (spatial) d dimensions being 
compactified, typically on a torus T^'^^ with radius R^ each way. This is called the bulk, as 
opposed to the brane on which the observable Universe lives. Since gravity experiments 
have not really probed the sub-millimetre 

The actual value of Newton's constant G^^'^'^^ in the bulk can be taken as large as 
0{1 TeV~^), but its value G^^^ in the effective 4-dimensional space at length scales ^ Rc 
is the extremely small one (9(10~^^) TeV~^ measured in gravity experiments. The two 
are related, using Gauss' Law M, by 



l2 



Kc 



1 2+d 



where M^pi^'^^ 



G^ 



{4+d) 
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2 + d 

denotes the Planck mass in the relevant number of dimen- 



sions. If M^pi~^^^ ~ 1 TeV, then i?, ~ io30/rf-i9 ^_g_^ d = 1, R^ ^ 10^^ m, which 
predicts deviations from Einstein gravity at solar system scales. Since no such effects 
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are seen, we are constrainedQ to take d > 2. For these values Rc < I mm and there is 
no conflict with known facts[^. The smallness of Newton's constant Gj^"* is thus a di- 
rect consequence of the compactification and hence there is no hierarchy problem in this 
theory^. 

In traditional Kaluza-Klein (KK) theories the mass-spectrum of non-zero KK 
modes is driven to the Planck scale Mpi . This problem is avoided in the AADD model 
by having the SM particles live on a 'wall' with negligible width (which we identify with the 

brane). While SM fields are confined to the brane, gravitons correspond to excitations 
of closed strings propagating in the bulk. The only new effects observable on the brane 
— to which corresponds the observable Universe — will be those due to exchange of 
gravitons between SM particles on the brane. 

To construct an effective theory in 4 dimensions, i.e. on the brane, gravity may be 
quantized taking the usual weak-field ('linearized') limit, assuming that the underlying 
string theory will ultimately take care of the well-known ultraviolet problems. Interactions 
of gravitons now follow from the (4 + (i)-dimensional Einstein equations in the compactifi- 
cation limit. Feynman rules to the lowest order in k = \J IGvrG^^ for this effective theory 



have been worked out in detail in Refs. and |]T0[- We make use of the prescriptions of 

Han, Lykken and Zhang in our work. For convenience, some of the relevant Feynman 

rules are listed in Appendix A. 

In this effective theory, the couplings of the gravitons to the SM particles will be 

suppressed by the Planck scale M^p^ ~ 1.2 x 10^^ GeV. This is offset, however, by the fact 

that, after compactification, the density of massive KK graviton states is very high, being 

r f4)i ^ (4) 

indeed, proportional to Mp( . The Mp( dependence cancels out, therefore, leaving 



an interaction of electroweak strength, whose scale is set by the bulk Planck's constant 
Ms ~ M^p^'^'' ~ 1 TeV, henceforth called the 'string' scale. A further assumption usually 



■^Of course, for d = 1, we have Rc oc Mg^: thus Ms > 10'' TeV would make d — 1 viable. However, 

then there are no interesting collider signals. 

'^A related problem, that of stabilization of the compactification scale, exists, however; this has been 

discussed in Ref. 0. In fact, this is the chief criticism of the AADD model. 
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made is that the brane itself is heavy and its vibration modes decouple from the processes 



under consideration. This corresponds to a static approximation for the brane 111 



The interactions generated by gravity between SM particles on the brane can be 
written mostly in terms of exchange of spin-0 and spin-2 gauge bosons. At the lowest 
order in /t, the spin-1 gauge bosons decouple [H, |n| from matter due to the diagonal 
structure of the energy-momentum tensor for matter fields on the brane. Interactions of 
the spin-0 dilaton (or 'radion', R) and the spin-2 gravitons(G') — or rather, their massive 
KK modes — have been described in Refs. ^ and |jlO[ and several phenomenological 
studies at the tree-level have been made using these interactions. The salient features of 
such studies are briefly surveyed in the next section. However, it suffices to note that 
almost all of these tend to provide lower bounds on the string scale M5. 

One-loop effects of the interactions of towers of KK states of the graviton and radion 
have not yet been investigated in detail. The pioneering work in this direction was in 
section 3.5 of Ref. |l^, in which a leading-order calculation of self-energy corrections 
to the mass of a scalar particle were worked out. A more detailed calculation of the 
anomalous magnetic moment of the muon was attempted in Ref. |]T2[. In the latter, 
the result turned out to be 'remarkably finite', though the actual numbers were rather 
disappointingly small. It was argued in Ref. that one-loop effects are worth calculating 
even though there may be other TeV-suppressed operators at the tree-level, simply because 
the one-loop effects are completely calculable, given the Feynman rules. One can add the 
argument that a study of one-loop effects brings out subtle features of the theory in a 
way that tree-level effects can never be expected to do. 

In the present work, therefore, we compute one-loop corrections to the masses of the 
electroweak gauge bosons W and Z, which are constrained by the famous p parameter |1^ . 
Extra contributions to the p parameter have been discussed in the context of AADD-type 



models in the literature |jTj], but not in the context of one-loop calculations. It turns 
out that the one-loop contributions from virtual radion and graviton states are strongly 
divergent in the ultraviolet (UV), a feature which is to be expected in an effective theory. 
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The use of the string scale Mg as a momentum cutoff has been suggested in Ref. 
(and elsewhere) and this leads to large contributions which actually grow with increasing 
Ms- Recognising that this is not a pathological feature but simply a manifestation of 
the fact that we have an effective (non-renormalisable) theory, we can therefore, use the 
experimental bounds on the p parameter to derive upper bounds on this scale M^, unlike 
all the tree-level processes considered hitherto, which are suppressed by powers of Ms- 

An interesting feature which emerges in our calculation is the presence of infrared 
(IR) divergences (for d = 2, 3 and 4) in the summation over masses of particles in the KK 
towers. These are due to the zero modes of the graviton (and radion) KK states and arise 
ultimately from the fact that in the IR limit when the graviton momenta are comparable to 
the mass gap between neighbouring excitations, the continuum approximation is invalid. 
One is therefore led to believe that the IR divergences are spurious and should cancel 
out of a complete calculation^ A useful suggestion made by the authors of Ref. |Ty] is to 
regulate these expressions by simply dropping the zero mode (whose coupling is Planck- 
scale suppressed) and start from the first massive KK mode. Cancellation of the IR 
divergence will therefore manifest itself as a cancellation of the dependence of the final 
result on the mass of this mode. We shall see that with the present level of knowledge, 
the presence of these IR divergences makes the p-parameter calculation non-predictive for 
(i = 2, 3 and 4. However, these same IR divergences allow us to make a concrete prediction 
for (i = 5 and 6, even without a complete knowledge of the relevant Feynman rules. 

The plan of this article is as follows. Section 2 consists of a micro-review of phe- 
nomenological studies of the AADD model carried out at the tree level. In the next 
section, we describe the basic formalism required to calculate the VF, Z self-energy correc- 
tions (and hence the p parameter) in the AADD model. Section 4 is devoted to numerical 
results from the calculation and their interpretation. It is in this section that we discuss 



^ One should also note that gravity is known to be infrared-safe, especially in the full (4+(i)-dimensional 
theory. Compactification corresponds to replacing certain components of momentum by masses and this 
cannot generate IR divergences. The authors are grateful to Ashoke Sen for bringing this fact to their 
notice. 
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the cancellation of IR divergences. We summarize our results in Section 5. Appendix A 
comprises of the relevant Feynman rules. The machinery of one-loop integrals using the 
ultraviolet (UV) cutoff Ms is described in detail in Appendix B. 



2. Phenomenological Studies 

Since gravitons (and radions) couple to any particle with a non-vanishing energy- 
momentum tensor, it is possible to make a variety of phenomenological studies of the 
interactions of their KK modes. Though these have not yet been completely explored, 
several important results are already available. These investigations can be classified 
into two types: those involving real KK graviton production, and those involving virtual 
graviton exchange. A real KK mode of the graviton will have individual interactions 
with matter suppressed by M^pj and will, therefore, escape any detector of terrestrial 
proportions. One can, therefore, see signals |]T5[ with large missing momentum and energy 
if an observable particle is produced in association with a tower of KK gravitons. Cross- 
sections for these depend explicitly on d, the number of extra dimensions. Each process 
can be used to obtain a bound on the string scale Ms for a given d. The most dramatic 



of these is Ms > 50 TeV for d = 2 and it comes from a study |T6[ of the flux of neutrinos 
from the supernova SN1987A. However, this astrophysical bound drops to about a TeV 
as soon as we go to d > 3. Most of the other processes lead to lower bounds of about 
1-1.1 TeV on the string scale for d = 2, and are even weaker for d > 3. 

Virtual (KK) graviton exchanges lead to extra contributions to processes involving 
SM particles in the final state and can be observed as deviations in the cross-sections and 
distributions. After summing over all the KK modes of the graviton, a tree-level Feynman 
amplitude involving graviton exchange is proportional to A/M|, where 

|A| ~ log(M|/s) for rf = 2 , (1) 
~ {d-2)-^ forrf>2. 
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to the leading order in s/Mg, where ^/s is the scale of the interaction T^|. It has been 
found convenient to absorb the magnitude of A into M^'*, so that bounds from virtual 
processes can be presented in a compact form. However, it should be borne in mind that 
the effective string scale Ms is, then, Ms = X^^^^Ms, and hence different from the actual 
Ms, such as is considered for real graviton emission. Once this slight abuse of notation is 



understood, each process can be used [0 to obtain a bound on Ms for a given sign of A. 
Among the most stringent of these bounds are Ms > 1 TeV for A = ±1 which come from 
studies of experimental data on dileptons and dijets |19| at the Fermilab Tevatron. 

A summary of the most stringent lower bounds from collider data on Ms, as well as 
the astrophysical bound from a study of the neutrino flux from the supernova SN1987A, 
is given in Table 1. 



Process 


Ms 


d^2 


d = 3 


(i = 4 


d=5 


d = 6 


pp e+e- 


- 1.0 


- 1.2 


1.0 


0.84 


0.76 


0.71 


pp dijets 


- 1.1 


- 1.3 


1.1 


0.92 


0.84 


0.78 


astrophysical bound 




- 50 


^ 4 


- 1 







Table 1. Illustrating lower bounds on the string scale Ms (in TeV) from dilepton jjq/ and dijet jj^/ data 
at the Fermilab Tevatron. The astrophysical bounds from a consideration J7^/ of the neutrino flux at the 
supernova SN1987A are also shown for purposes of comparison. 



4. Calculation of Self-Energies 

In electroweak theory, the p parameter is defined by [jl3 



Ml cos2 Ow 



(2) 



to all orders in perturbation theory. The most convenient formalism for the study of this 
parameter is that given by Peskin and Takeuchi who define p—1 = aT, where a is the 
fine structure constant and p is calculated in the zero momentum limit. The T-parameter 
defined above is one of the so-called oblique parameters and is zero in the Standard Model 
at tree-level. Non-vanishing values of T measure, therefore, one-loop corrections to the 
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W, Z masses in the SM and/or any new physics effects in these masses. The current 
experimental bounds on the (dimensionless) T-parameter are T = —0.21 ± 0.16 (+0.10). 
In terms of the self-energy corrections n(p^) of the gauge bosons and Z^, the 



(3) 



T-parameter is pO[] 

Ati r 

'^ = J7r ^ww{0)-Uzz{0)cos^ew 

where p is the propagator momentum and the self-energies are calculated in the limit 
0. This will have small contributions in the Standard Model. The excess contribution to 
the T-parameter in the AADD model (henceforth denoted simply 6T) can be calculated 
simply by evaluating the self-energy corrections of the W and Z bosons due to graviton 
(radion) loops. 

(a) 



W ^ W* 




w 



(b) 



f,W,Z,H 



W 



K 



w 



(c) 



w 




G*/R* 



i G*/R* 



W 




V 2 

w ^ w 

Figure 1. Feynman diagrams corresponding to self-energy corrections of the W-boson in the AADD 
model at the one-loop level. G„ and _R„ correspond to graviton excitations with spin 2 and respectively. 
There will be a similar set of diagrams for the Z-hoson. 

The Feynman diagrams corresponding to the self-energy of the VT-boson are shown in 
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Figure l{a-d). Since the low energy theory is just hnearized Einstein gravity, we expect 
it to be non-renormahzable. Hence, rather than work with a whole family of effective 
operators Q, we work with the unrenormalized vertices as given in Ref. |[10| and it is 
therefore necessary to include tadpole graphs such as those of Figure 1(6) and 1(c). These 
can have all kinds of particles in the loops, including (6) all SM particles and (c) the 



quanta of linearised gravity. Moreover (as pointed out in Ref. [10|) it is necessary to 
include 'seaguU'-type diagrams such as that in Figure l{d). It is important to note that 
while the diagrams in Figure l(a-c) involve 0{k,) vertices, the diagrams in Figure l{d) 
involve O(k^) vertices. Unfortunately, Feynman rules are not readily available in the 
AADD model for (i) the self-couplings in the pure gravity sector, which are required for 
the evaluation of Figure 1(c) and {ii) the full set of 0{k,'^) vertices. Thus, any calculation 
of the self-energies at this stage must be incomplete to this extent. However, as we shall 
show, it is still possible to obtain some meaningful results, and this is the theme of this 
work. 

Evaluation of the diagrams in Figure 1(a) and (6), using the Feynman rules given 
in Appendix A, is a long and tedious process. Some of the formalism is developed in 
Appendix B. The final result for H^|^/^^(0) takes the form 



H^^^^(O) ^ ^ t dx x-^^^l^ 



(4) 



where x = M^/Mg. The functions J^'^(x), which arise from evaluation of the loop 
integrals in the case of graviton (G) and radion [R) respectively, are defined in Appendix 
B. In the above formula xq is an IR cutoff which is taken to be xq = {RcMs)~'^ following 



the suggestion of Ref. [TQI- It corresponds to the mass of the first massive KK state of 
the graviton/radion. Using the relation between Ms and rf, it follows that 

which is clearly minute for Ms ~ 1 TeV and approaches unity when Ms — > M^p} . 



^such as, for example f fWW-ty^e and f f Z Z-type terms; there are many other possibiUties. 
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(6) 



Similarly, the final result for 11^^^^ (0) takes the form 

^ ' 7207r (2v/^)'^ r(rf/2)ixo 
where, as before, the functions Jz'^i^) are defined in Appendix B. Substitution of Equa- 
tions and in Equation now yields the value of 6T for a given d and value of M5. 
Since the functions Jw,zi^) extremely complicated (see Appendix B), the integrals in 
the above formulae are evaluated numerically. 

Evaluation of the Feynman diagrams in Figure 1(c) and {d) is not directly possible 
because the Feynman rules are still not available. However, we can make an approximate 
estimate of the 'seagull'-type graph using the Feynman rule given in Eqn. (84) of Ref. 
p!0| for a pair of (spin-2) gravitons coupling to a pair of scalars. Assuming that the 
longitudinal components of the W and Z bosons are given by the Goldstone equivalence 
theorem, we can evaluate the 'seagull' graph. This adds to the function J^^(x) by a term 
which we denote K^^i^)- The full form of K^^i^) is given in Appendix B. We then 
make the assumption that the two contributions due to the transverse modes are each 
equal to that due to the longitudinal mode. This corresponds to multiplying Ky^ ^{x) by 
a factor of 3. We next parametrize the error due to these approximations by introducing 
an unknown parameter ^, so that the extra contribution reads 

JwA^) JwA^) + 3 e K^,zi^) (7) 

We next note that the diagrams in Figure 1(c) have the same particles running in the 
loop as the corresponding ones in Figure l{d), and therefore, should be proportional to 
the same loop integrals. In fact, it will not matter which of the graviton states contributes 
to the loop integral, since the propagators are same, each having an identical tower of 
states. Of course, there will be explicit dependence on the graviton/radion in the vertex 
and in the propagator connecting the loop with the WW vertex. However, to a first 
approximation, we assume that all these can be lumped into a common contribution, also 
proportional to K^^^i^x). All the extra contributions due to the diagrams in Figure 1(c) 
are thus absorbed into the unknown parameter ^. Any momentum dependence arising 



from the vertices and propagators contributing to ^ can be removed by defining ^ to be 
an average value. We therefore, do not add any extra contribution to Jj^ 



While the above procedure does not sound very rigorous, it is the best that can be 
done until all the relevant Feynman rules become available. Calculation of these is a non- 



trivial task and it may be some time before they are available in the literature [^. In 
principle, the value of ^, which parametrizes our present state of ignorance of the effective 
theory to can depend explicitly on Ms and on d. However, as we show in the 

next section, any such dependence is very weak, and therefore, we obtain an a posteriori 
justification for the approximation made here. 

It is clear from the above that we have little or no a priori knowledge of what the 
parameter ^ should be, except that it should be of order unity (which follows from di- 
mensional arguments). However, it is here that the IR divergences come to our help. We 
shall see in the following section that gravitonic contributions to bT are IR divergent for 
= 2, 3 and 4. It has been argued above that these IR divergences are spurious and hence 
we demand that they should cancel out when we sum over all diagrams. Since Jwzi^) 
and K^^^{x) are individually IR divergent, we can tune the value of ^ to get cancellation 
of these divergences. We then claim that this is the value of ^ which would be obtained 
from a proper knowledge of the Feynman rules in this modelQ. The exact procedure is 
described in the following section. 

4. Numerical Results and Discussion 

Numerical evaluation of the integrals in Equations (4) and (6) is not completely 
straightforward. To see this, consider, for example, Ms = 1 TeV and d = 2, in which 
case, the IR cutoff comes out to be xq ~ 10^'^^. The integral runs, therefore, over several 
orders of magnitude and, being IR divergent, receives its principal contribution from the 
smallest values of Xq. It is, therefore, convenient to smoothen the integrand by the trans- 
formation ^ = — logx, which allows the integration to be done using standard numerical 
^This is in the same spirit as the determination of, for example, the three-gluon vertex in QCD from 
considerations of gauge invariance 
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techniques. For d < 4, the IR cutoff yo = — logxo grows with Ms (corresponding to the 
increase in energy gap between the two lowest KK states), so that, for very large values of 
Ms ~ M^p! (which means xq — 1) the range of integration gets pinched off. However, the 
pinch-off is very sharp and only takes place when xq — 1 to great precision. Of course, 
such large values of Ms are uninteresting from the experimental point of view, but it is 
reassuring to note that even if we get unacceptably large values for the p-parameter at 
low values of Ms, there is always the possibility of having extra dimensions compactified 
to the Planck length. 

In their leading order calculation of one-loop effects in the AADD model, the authors 
of Ref. |]TU[ have shown that these corrections to the scalar mass can be written in terms 
of (dimensionless) loop integrals, of which two are 

h = t dx x-^+^/^ , h= C dx ^ (g) 

J XQ J Xo 

where Xq is the IR cutoff described above. On evaluation these come out to be 

Now, these are clearly IR divergent for d < 2 and d < 4 respectively. Even the use of the 
IR cutoff Xo leads to extremely large values of these integrals, since it is clear that xq is 
an extremely small quantity (unless Ms — >■ Mp^ . 

We now adopt the philosophy that the infrared divergences are spurious and must 
cancel out of a complete calculation. The main motivation for this — explained above — 
is that gravity is well-known to be infrared- safe, in fact, more so in dimensions greater 
than four, and this feature is expected to be preserved under compactification. The 
origin of the IR divergences may be attributed to a breakdown of the mass-continuum 
approximation in the low- momentum limit. Since we cannot make a complete analytic 
calculation and check the cancellation explicitly, the best we can do is to tune the value 
of the unknown parameter ^ and see if we can achieve cancellations in the different cases 
for d and Ms- Our results are shown in Figure 2, for Ms = 1 TeV, Mh = 250 GeV and 
d = 2,3 and 4. 
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Figure 2. Illustrating the cancellation of IR divergences in 6T for d = 2,3 and 4- Note that these 
occur for the same value o/^ in each case. This is also insensitive to variations in Ms and Mh- 

In Figure 2, it may be seen that a sharp cancellation occurs for a specific value 
^ ~ 3.2836. Amazingly enough, we get the same value for all the three values of d. 
We also find that the value of ^ for which the cancellation occurs is almost completely 
insensitive to the value of Ms, while it has a very weak dependence on the mass of the 
Higgs boson (0.2% as Mh varies from 100 GeV to 700 GeV). This constant value of ^ may 
be attributed to the fact that it is almost entirely made up of constant or averaged factors 
in the unknown coupling, including k, Rc and the gauge boson masses, apart from purely 
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numerical factors. We believe that a complete calculation, when available, will predict 
precisely this value ^ ~ 3.2836. It is reassuring to note that the value of ^ is indeed of 
order unity as we expect. 

Obviously, when such fine cancellations take place, it is not possible to predict the 
value of 5T for o? = 2, 3 and 4, since any result can be obtained by making small changes 
in ^. However, in view of the fact that ^ is almost constant for o? = 2, 3 and 4, we feel 
emboldened to extrapolate this value to the (IR-finite) cases d — b and 6 as well. We thus 
feed in this value of ^ for the cases d — h and 6, and are able to make concrete predictions 
for 5T in the AADD model. Figure 3 shows our results for 5T as a function of Ms (a) in 
the range of interest for collider experiments and (6) in the region where the UV-divergent 
integral gets pinched-off as ^ 1. As before, we set Mh — 250 GeV, though there is 
little variation of the results with this parameter. 




0.4 0.8 1.2 1.6 2 0.4 0.8 1.2 1.6 2.0 

Ms (TeV) Ms (10^^ TeV) 

Figure 3. Variation in ST in the AADD model with increasing Ms, setting ^ = 3.2836. In (a), the 
dashed lines correspond (from bottom upwards) to the experimental upper limits at la, 2a, 3a respectively. 
The pinching-off of the UV-divergent integral at the Planck scale is illustrated in (b). 

As Figure 3 shows, the gravitational contribution to the T-parameter shows a non- 
decoupling behaviour, though, of course, it finally falls to zero when Ms M^ . The 
non-decoupling behaviour, we feel, is an indication that the theory is an effective one. In 
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fact, for d = 5, it appears that any value of Ms in the range displayed in Figure 3(a) is 
ruled out by experiment — which forces us to have Ms ~ Mpi , as in Figure 3(6). For 
d = 6 we obtain an upper bound of about 725 GeV at 95% C.L., which is barely allowed 
by experimental data (see Table 1). A slight improvement in the data is likely to rule this 
out and this would force Ms ~ Mpi for d = 6 as well. 

The actual scenario need not be so gloomy, however, since the effective theory of 
Refs. and is only an approximation and the underlying (string?) theory is believed 



to be finite. At the present juncture, it is not possible to make calculations in the exact 
theory, and it would, therefore, be necessary to construct a phenomenological extension 
of the effective theory, by adding extra (non-renormalizable) operators with unknown 
coefficients^ to the Lagrangian The coefficients may then be adjusted to produce 

acceptable values of the p (and T) parameter, even for d = 5 and 6. This is a messy 
business which smacks of fine-tuning and will not be attempted in this work. 

Finally, it is worth noting that in all these calculations, the contribution due to the 
radion, which has much lower levels of divergence, both in the UV and IR regions, turns 
out to be negligible in comparison to that of the graviton. Its inclusion in all the above 
formulae is more for the sake of completeness that for numerical accuracy. 

6. Conclusions 

In this work, we have performed a careful study of the W and Z boson self-energy 
corrections using the effective theory of KK graviton (and radion) interactions developed 
in Ref. |]10[ from the original suggestion of Refs. JI], ^. Since these are one-loop calcu- 



lations, the results have several interesting features, which do not show up in tree-level 
calculations. The results turn out to be strongly divergent, both in the UV and IR 
regimes. Using a cutoff at either end, as suggested in Ref. can lead to large values 



of the T-parameter (and hence of the p parameter). Knowing gravity to be IR-safe, we 
demand cancellation of the IR divergences and use this to determine the unknown which 



^This would be analogous to the twelve 0{p^) operators in chiral perturbation theory. 
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parametrizes our ignorance of the full set of Feynman rules in this model. Interestingly, 
we find that ^ is independent of Ms and the number of extra dimensions, and hence can 
be used to make concrete predictions for 5 and 6 extra dimensions. We then find that for 
(i = 5, Ms is driven by the T-parameter constraint to the Planck scale, while for d = 6, 
a small window in Ms around 725 - 750 GeV is still viable. We are unable to make any 
predictions for < 5. 

We wish to emphasize that while our results do seem to show that models with large 
extra dimensions do not work well for c? > 4, all that this means is that the formalism of 
Ref. |]10| does not work well for > 4. A deeper (perhaps finite) theory should certainly 



provide better insights and more acceptable values of the p-parameter. It might also be 
possible to build a phenomenological theory by adding non-renormalizable operators to 
the effective Lagrangian, but that lies outside the scope of this work. 
Note added: After this work was completed, we received the preprint in Ref. where a 
similar calculation has been done, using a purely analytic approach. The rather striking 
difference in our results from theirs is probably due to the fact that tadpole diagrams, 
as in Fig. 1 (6) and (c) have not been taken into account in Ref. |]2S|. We thus believe 



that the calculation in Ref. p5|, though elegant, is incomplete and should not be taken 



as evidence for a decoupling behaviour in Ms-, at least for the p parameter. 
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Appendix A: Feynman rules 



The Feynman rules relevant for the calculation of the p parameter are as follows. 
The vertex for a pair of gauge bosons coupling to the graviton are given by an effective 
Lagrangian 



The graviton coupling Tq^'^{p, q) is given by 



where n = ^/IQnGjq and 



(A.l) 



(A.2) 



(A.3) 



(A.4) 



In the above equation, both W/Z momenta are assumed to be directed towards the vertex. 
The radion coupling is given by 



rppa 



I 3(2 + d) 

The propagator for a graviton with momentum k is 

iPpypai}^^ 



fc2 - M2 + ie 



where the polarisation sum P^vpa is given by 



pupa 



1 / kfj^kp 
+ ^(Vpp 



Vpu 



M2 

n 

kpkcr 

M2 

n 

kpki^ 



Vup 
Vpa 



kiikfj 

n 

kj/kp 

n 

kpkcr 



(A.5) 



(A.6) 



(A.7) 
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Finally, the propagator for a radion is given by 

i{d-l) 



(A. 



A;2 - M2 + ie ' 

where the extra factor [d — 1) arises from the sum over different modes of the radion in 
each dimension. 



In this Appendix we have closely followed the notation and conventions of Ref. |T0 

Appendix B: Loop Integrals 

The momentum integrals arising in one-loop calculations involving KK graviton and 
radion modes are calculated using an ultra-violet cutoff Ms, where Ms is the 'string' scale. 
We use the well-known notation of 'tHooft and Veltman and Passarino and Veltman [53 



to describe momentum integrals of products of propagator functions. Then, the scalar 
one-point function turns out to be 

Aim') = ^ 



where 



TT^ k'^ + rri' + ie 
Ml a{^) (B.1) 



A{x) = l-x log(l + -) (B.2) 

X 



As in Refs. E^], all integrals are defined in Euclidean space. 



For the calculation of T, we require to calculate the boson propagators in the zero 
momentum limit. Thus, for the scalar two-point function, we define 

d^k 1 

p-^O J TT^ [fc2 _|_ iTjSj |-^^ _|_ py _|_ 



B„(...x.) = B.K.,„^;0) = lim/ 



-1 



(Xi - X2) 

where xi^2 = ^12/^3- then in a position to define 



A{xi) - A{x2) 



(B.3) 



d 

BQ{ml,ml]0) = lim — 5o(m?,m^;p2) 
P^^o op 

= -l^B'^(xuX2) (B.4) 
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which turns out, on evaluation, to be 

^o(a;i,a;2) = ^ 

XiX2{Xi - 



(B.5) 



X 



Xi 



Xi- X2- X1X2 l0g( ) - (1 + X1X2) {xiA{x2) - X2A{xi)} 

X2 



The vector two-point function is defined to be 



B^{m^,m2]p) = J ^ 



[k^ + mj] [{k + py + ml] 
Pi, Bi{mj,ml;p'^) , 



(B.6) 



which leads to 



Bi{xi,X2) = Bi{ml,ml;0) 
_ 1 
~ ~ 2 



Bo{xi,X2) + {Xi - X2) Bq{xi, X2) 



The tensor two-point function is 

d*k 

B^^{m{,ml]p) 



k^ kjy 



TT^ [A;^ + mf ] [{k + + ml] 
Pi,Pu B2i{ml,ml;p^) + 5^^ B22{ml,ml;p^) , 



(B.7) 



(B.8) 



and we define 



B22{ml,ml-Q) = Ml B22{xi,X2) . 
This function can be evaluated as 



B22{Xi,X2) = - 



A{xi) - 2X2 -Bo(Xi,X2) + {xi - X2) Bi{xi,X2) 



(B.9) 



(B.IO) 



A similar expression can be derived for B21, but it is not relevant for the present calcula- 
tion. 

Now, in terms of these functions, we can write the integrands relevant for 111^1^(0) as 



J^{x) = Ig{Xw,x 



(B.ll) 



x 



283 - 12 XwA{Xw) - 6 XzA{Xz) - 4 XhA{Xh) 



+ 8j2Cf XfA{Xf) 
f 
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where — Mi^ /Mg and the colour factor C/ = 3 for quarks and 1 for leptons. Similarly 



w 



16 Id{Xw,x) 
6 



(B.12) 



X 



295 + 24 XwA{Xw) + 6 XzA{Xz) + 2 XhA{Xh) 



+ 8j2Cf XfA{Xf) 
f 



Similarly, we can write the integrands relevant for Hzzi^) as 



J|(x) = -Ig{Xw,x) 



(B.13) 



SXz 
2x 



283 - 12 XwA{Xw) - 6 XzA(Xz) - 4 



+ 8 ^ 
/ 



AX, 



4 /o(Xiy,a;) 

--(295 + 24 XwA{Xw) + 6 ^^^(Xz) + 2 



(B.14) 



+ 8 XfA{Xf) 

f 

The /-functions, which arise from evaluation of the diagram in Fig. 1(a), can be 
written as 



Id{xi,X2) = ^ A{xi) + {3Xi - X2) Bo{xi,X2) + ^{Xi - X2) Bi{xi,X2) 



(B.15) 



for the loop integral involving radion exchange, and 

fl{xi,X2) + f2{xi,X2) A{xi) + f3{xi,X2) A{X2) 

+ Mxi,X2) 522(a;i,a;2) 



Ig{xi,X2) = \ 



(B.16) 

for the loop integral involving graviton exchange. The functions fi{xi,X2) are given by 



fl{xi,X2) = {X1-X2) 



12x1 - QOxixl + 4(24a;? - 19a;i - ll)a;2 



-{A8xl + 98x1 + 64xi + 33) 



(B.17) 
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f2{Xi,X2) 



48x1X2 12x2 + X1X2 — 3x1 



(B.18) 



/3(a;i,X2) 



12 xl - 47x1X2 + 10x1X2 - 4xi 



(B.19) 



/4(a;i,X2) 



48xi(xi— X2) X2 + 8x1X2 + Xl 



(B.20) 



Finally, the 'evaluation' of seagull diagrams leads to the functions 



28.9 \2 + 21x + 5Ax^A{x) 



X 




72.0 — [2 + 21x + 5Ax^A{x) 



(B.21) 



X L 
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